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Adiabatic time-dependent density functional theory fails for excitations of a heteroatomic molecule 
composed of two open-shell fragments at large separation. Strong frequency-dependence of the 
exchange-correlation kernel is necessary for both local and charge-transfer excitations. The root 
of this is static correlation created by the step in the exact Kohn-Sham ground-state potential 
between the two fragments. An approximate non-empirical kernel is derived for excited molecular 
dissociation curves at large separation. Our result is also relevant for the usual local and semi-local 
approximations for the ground-state potential, as static correlation there arises from the coalescence 
of the highest occupied and lowest unoccupied orbital energies as the molecule dissociates. 



PACS numbers; 31.15.Ew,31.10.+z,31.25.Jf 

I. INTRODUCTION AND MOTIVATION 

Time-dependent density functional theory 
(TDDFT) has seen a surge of applications in 
recent years [E S IS S IE B|j with calculations of 
electronic excitations in a wide variety of systems of 
physical, chemical, and biochemical interest. Errors 
in TDDFT excitation energies are typically within 
a few tenths of an eV, while implementations can 
scale with size even more favorably than that of 
time-dependent Hartree-Fock. In principle, TDDFT 
yields exact excitation spectra: the errors arc due 
to approximations to the exchange-correlation (xc) 
functionals. Approximations enter at two stages 
in usual linear response calculations: first, in 
the xc contribution to the one-body ground-state 
Kohn-Sham (KS) potential out of which bare exci- 
tations are calculated. For example, its too rapid 
asymptotic decay in LDA/GGA causes problems 
for high- lying bound states 0, H, S Second, 
the dynamic xc kernel must be approximated; this 
corrects the bare KS transitions towards the true 
transitions. Almost all calculations employ an 
adiabatic approximation, which entirely neglects 
frequency-dependence in the kernel. Its remarkable 
success for many excitations is not well understood. 
Understanding where it fails is important in con- 
structing more accurate functional approximations. 
The lack of frequency-dependence is why adiabatic 
TDDFT (ATDDFTWails for states of double 
excitation character [iTllTHlllTl Hsf . 

Molecular dynamics simulations, atom trapping, 
and photoassociation spectroscopy all need accu- 
rate calculations of excited molecular dissociation 
curves El IH El E HI 113, but these are 
challenging for approximate methods. Even ground- 
state dissociation is difficult Closely linked 
are the difficulties in obtaining long-range charge- 
transfer (CT) excitations. The latter are notoriously 



underestimated in ATDDFT IH HI H IH HI IH 
I2I I2I Is^, "sTl , yet are vital in biological and chem- 
ical systems large enough that density functional 
methods are the only feasible approach. 

We show here that a strong frequency-dependence 
in the kernel is needed near every excitation of a 
molecule composed of two different open-shell frag- 
ments at large separation: The ubiquitous adiabatic 
approximation fails for all excitations. We show this 
is due to static correlation in the KS ground-state, 
leading to every single excitation being almost de- 
generate with a double excitation. We resolve what 
the frequency-dependence of the exact xc kernel is 
for the tails of excited dissociation curves, and derive 
an approximation based on this. Our kernel yields 
both charge-transfer and local excitations, and con- 
tains the correct dissociation limit and leading-order 
polarization. It bootstraps an adiabatic approxima- 
tion, and yields results to the same good level of 
accuracy that ATDDFT has for more common local 
excitations: that is, our kernel undoes the problem- 
atic static correlation in the KS system. 

A crucial role is played by the step structure in the 
exact ground-state KS potential that develops as a 
heteroatomic molecule dissociates |3l IH |3 Isl 
Is^ . The step has size equal to the difference in the 
ionization potentials of the two fragments, so aligns 
the highest occupied molecular orbitals (HOMOs) 
of each. Its consequences for TDDFT have only be- 
gun to be explored [23 . Is^ . Recently we showed 
that near the lowest CT excitation on each fragment, 
the step imposes a st rong frequency-dependence on 
the exact kernel |3j l38| . Here, we go much fur- 
ther, showing the dramatic frequency-dependence 
required for both higher CT and all local excitations. 
Local and CT excitations are intimately entangled 
by the step in the KS system and we construct the 
kernel that disentangles them. 

It is important to note that although local and 
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semi-local approximations to the ground-state po- 
tential (ie LDA and GGAs) lack the step, static 
correlation still features, since in these cases, the 
HOMO and lowest unoccupied molecular orbital 
(LUMO) become degenerate as the molecule disso- 
ciates. So, the construction of our xc kernel has 
relevant implications for the case of the usual ap- 
proximate ground-state KS potentials too. 

In Sec. ^ we briefly review linear response the- 
ory in TDDFT. We discuss its apphcation to CT 
states of a long-range molecule and point out a fun- 
damental difference in the analysis when the closed- 
shell molecule consists of open-shell fragments rather 
than closed-shell fragments. We summarize the re- 
cent result of Ref. on the structure of the exact 
kernel that captures the lowest CT excitations ex- 
actly in the case of open-shell fragments in the large- 
separation limit. In Sec. IIIII we extend the analysis 
to all excitations of the molecule, showing how both 
local and CT excitations require strong frequency- 
dependence in the xc kernel that is lacking in all 
approximations to date. We derive an approxima- 
tion to the kernel and demonstrate it on a simple 
molecule. Finally, Sec. IIVI contains a discussion on 
implications of our result. 



II. TDDFT LINEAR RESPONSE AND 
LOWEST CHARGE-TRANSFER STATES 

TDDFT linear response 0, 13 proceeds by cor- 
recting the KS single excitations towards the true 
excitations, (which may be mixtures of any num- 
ber of excitations) through the Hartree-exchange- 
correlation kernel /Hxc[«o](r, r', t^^) = l/|r — r'| + 
fxc[no]{'ic , r' , uj) , a functional of the ground-state 
density no(r). The xc kernel is the frequency- 
time Fourier-transform of the functional deriva- 
tive of the xc potential, /xcij^oKr, r', t — t') = 
dvxc{j^t)/Sn{r't)\no- The density-density response 
function of the KS system Xs smd that of the in- 
teracting system x are related through: 

r'M = Xs~'H-/HxcH. (1) 

Usually a matrix version Q is used to obtain ex- 
citation energies and oscillator strengths. The vast 
majority of calculations utilize an adiabatic approx- 
imation for the xc kernel, i.e. one that is frequency- 
independent. Typically the results are accurate to 
within a few tenths of an eV (although there are 
exceptions, see eg. H, [23 ) ■ 

When excitations are well-separated, the matrix 
may be simplified in a "single pole approximation 
(SPA)" m H 113 that gives the true frequency w 
as a shift from a KS transition: for a closed-shell 



molecule, 

UJ = UJq + 2[q\fHxc{'^q)\q] ■ (2) 

Here, q = (i, a) is a double-index labeling the sin- 
gle KS excitation from an occupied orbital 4>i to 
unoccupied 0^, ujg = Ca - and [q\fnxci(^)\q] = 
J drdr'0*(r)0a(r)/Hxc(r, r', c.)</.,(r')<^:(r'). 

The SPA is also an important tool in understand- 
ing the TDDFT corrections to the bare KS ener- 
gies. Applying Eq. Q to the problem of lowest 
charge-transfer states in a long-range molecule p5| . 
shows immediately that the correction from the ker- 
nel, [q\f-a:Kci'^q)\q], vanishes, because there is van- 
ishing overlap between the occupied orbital (r) on 
the donor and the unoccupied 0a (r) on the accep- 
tor, in the limit of their large separation R. The 
TDDFT value of Eq. ^ then reduces to the bare 
KS eigenvalue difference, i.e. the acceptor LUMO 
orbital energy minus the donor HOMO orbital en- 
ergy. The latter is equal to minus the ionization 
energy of the donor, In, while the former is the KS 
electron affinityHllHllillll, Aa,s = Aa- Aa,xc, 
not the true electron affinity of the acceptor, Aa- 
The resulting frequency, uj = Ijj — Aa,s, is a se- 
vere underestimate to the true CT energy: aside 
from the usual local/semi-local approximations un- 
derestimating the ionization energy, it lacks the xc 
contribution to the electron affinity, as well as the 
tail 25, The exact CT energy is 

Id — Aa — 1/R in the large-R limit. 

When the (closed-shell) long-range molecule con- 
sists of two open-shell fragments |23 (labelled as 
1 and 2 below), this argument needs to be revis- 
ited because neither the HOMO nor LUMO are lo- 
calized on one of the two fragments. In the large 
separation limit, the exact KS potential contains a 
step [33, IS 42] that exactly aligns the individual 
homo's of the fragments, rendering the HOMO 
of the long-range molecule, (/)o, to be the bond- 
ing molecular orbital composed of the HOMO's of 
the individual fragments, 0o — (0ff,i + 0h,2)/V2- 
The LUMO 0g is the anti-bonding molecular or- 
bital. (See Fig.lR and also e.g. Fig. 1 of Ref. J^], 
Fig. 7 of Ref. Ha, and Fig. 1 of Ref. 0). The 
lowest CT states lie in the subspace formed by the 
HOMO and LUMO, as shown by the explicit diag- 
onalization of the Hamiltonian in Ref. [2J|. Turn- 
ing to TDDFT, Eq. @, the bare KS energy differ- 
ence, ujq, vanishes as an exponential function 
of the separation of the two fragments, i?, since 
the antibonding LUMO and bonding HOMO are 
separated in energy only by the tunnel splitting 
through the barrier created by the step. The cor- 
rection from the xc kernel in Eq. |21 is no longer 
zero, as these orbitals have significant overlap in 
both the atomic regions. However any adiabatic ap- 
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proximation to the kernel will lead to drastically in- 
correct CT energies: The underlying reason is be- 
cause the adiabatic approximation entirely neglects 
double-excitations [llLll^ . which are also absent in 
the non- interacting response function Xs- Yet the 
double excitation to the antibonding orbital is criti- 
cal in this case to capture the correct nature of the 
true excitations in this subspace p3| (as well as the 
Heitler-London ground-state). Refs. [lllIT^ showed 
how the mixing of a double excitation with a single 
can be incorporated in TDDFT via a dressed SPA, 
uj ~ LUq + 2[(7|/hxc(w)|(;]: here the kernel has strong 
frequency-dependence. In Ref. the form of the 
exact non-adiabatic xc kernel was uncovered in the 
limit of large separation and in the limit that the 
CT excitations have negligible coupling to all other 
excitations of the system. This was obtained by ex- 
amining the structure of the response functions x 
and Xs in this subspace and comparing with the ex- 
plicit diagonalization of the Hamiltonian. The result 
was: 

oj[q\fH^c{oj)\q] = S + + — 3 

4 uj^ — UO1LO2 

for the KS transition q =bonding HOMO— >anti- 
bonding LUMO orbital. Here, wi = I2 — Ai — 
1/R and uj2 ^ h - A2 - l/R, where A = A + 
A^^'°^, with the xc contribution to the electron 
affinity of each fragment approximated as ^^c'""' = 
-Jd^rJ d3r>i/(r)20^(r') V|r - r'|, and (j)H is the 
HOMO of the appropriate fragment. Finally, S — 
{lui - W2)/2, and w - e"^ is the HOMO-LUMO 
sphtting (the bare KS frequency) [Mill- Notice the 
strong non-adiabaticity in the exact kernel, manifest 
by the pole in the denominator of Eq. Q . 

The strong-frequency dependence in this expres- 
sion arises due to static correlation in the system: 
the near-degeneracy of the antibonding orbital leads 
to the KS ground-state Slater determinant being 
near-degenerate to two others. This feature makes 
heteroatomic dissociation in the ground and lowest 
CT states resemble homoatomic dissociation in some 
respects ; problems associated with ground-state 
homoatomic dissociation are now well-known |43|. 
In the present paper we explore the consequences 
of the static correlation on all other excitations of 
the heteroatomic system, and model the required 
frequency-dependence in the kernel. Every single ex- 
citation out of the occupied bonding orbital is almost 
degenerate with a double excitation, where another 
electron is excited from the bonding orbital to the 
antibonding orbital . This near-degeneracy leads to 
strong-frequency dependence in the xc kernel rele- 
vant to all excitations of the molecule. 

It is very interesting to note that although 
the usual local/semi- local approximations 



(LDA/GGA's) to the ground-state potential 
do not contain the step, the HOMO and LUMO 
are nevertheless also delocalized over the long-range 
molecule. We shall return in Sec. II VI to a discussion 
of this. 



III. LOCAL AND HIGHER CHARGE 
TRANSFER EXCITATIONS 

Consider first a model molecule composed of two 
different "one-electron atoms" . At large separations, 
the true Heitler-London ground-state (Fig.QjV), has 
one electron on each atom (in orbitals (jji and 4>2 re- 
spectively). On the other hand, due to the alignment 
of the atomic levels caused by the step in the KS po- 
tential (Fig.^B), the KS ground-state is the doubly- 
occupied bonding orbital, 00 = (0i + 02)/\/2, with 
each electron evenly spread over each atom. Gener- 
ically, higher excitations of the atoms do not coin- 
cide; for our example, atom 2 has one higher bound 
excitation, whereas atom 1 has none. 




FIG. 1: A: Model of long-range molecule. 
The Heitler-London singlet ground-state is 
((^i(r)(?!)2(r')-H<;/'2(r)<^i(r'))/v^. The dashed curve 
is an excitation on "atom" 2. B: KS potential, (which 
eventually returns to zero on the right). The KS 
ground-state is the doubly-occupied bonding orbital (jio- 

The vertical arrow in Fig. illustrates a local 
excitation on atom 2, (/)2 4>2* ■ If we attempt to 
describe this in the KS system, we immediately see 
a problem: Excitations must occur out of the oc- 
cupied bonding orbital 0o (Fig. QjB), so excitation 
leaves the molecule with "half an electron" on atom 
1 and "one and a half electrons" on atom 2. This sin- 
gle excitation, 0o ~^ 4'2* is a poor description of the 
excitation of the true system. The true excitation is 
in fact a linear combination of this single excitation 
and a double excitation, (^Oj'/'o) ~^ {4>2*t4>q) where 
the other electron occupying the bonding orbital is 
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excited to the antibonding orbital (po . (We assume 
for now the other excitations of the system are far- 
ther away in energy). Because the transition fre- 
quency to the antibonding orbital is the tunnel fre- 
quency between the atoms, it is exponentially small 
as a function of their separation. Thus, every single 
excitation of the system ((po (pa) is almost degen- 
erate with the double excitation (0Oj (po) ^ i'Pa, 4>a)- 
This feature is a signature of static correlation in the 
KS ground state. Since double excitations are miss- 
ing in the KS linear response it is the job of the 
xc kernel to fold them in (see shortly). 

It is instructive first to diagonalize the interacting 
Hamiltonian H in the 2x2 KS basis, whose KS 
energies differ only by the tunnel splitting(~ e^^): 

1>5 = ((/.o(r)02-(r') + 02*(r)0o(r'))/\/2 
$D = (Mr)'/'2-(r') + 02*(r)0^(r'))/V2 (4) 

H sums the kinetic, external potential, and electron 
interaction (Ke) operators, and rotates the pair into: 

*2- = (0i(r)02.(r')+02.(r)0i(r'))/\/2 
*CT2. = (02(r)<^2-(r') + '/'2-(r)02(r'))/V2 (5) 

Within the truncated space, Eqs. © are eigenstates 
of the interacting system: ^'2* is a local excitation 
on atom 2, while ^'ct2* is an excited CT from atom 
1 to an excited state of atom 2 (dotted arrow in 
Fig. n^). Because these states are paired together 
due to the bonding orbital nature of the KS ground 
state, the two qualititatively different excitations 
appear together in the exact TDDFT, as we will 
see shortly. The eigenvalues of the diagonalization 
give their approximated frequencies, after subtract- 
ing the expectation value of H in the Heitler-London 
ground-state, Ehl-, (see also Refs. [TTll2^): 

<^2* = £2* ~ £2 

0OCT2' = e2-~ei-4'*^-l/i?, (6) 
to leading order in R, where 

A^P = - j j dr'(|02.(r)n(/-2(r')P 

+ {r)<Pl{v')c^2' (r')02(r)) 14c(r - r')(7) 

This approximates the xc part of the excited electron 
affinity A*^^*^: generally for an iV-electron species, 

A* ^Ea{N)~E*{N + l)=Al+A*^^ (8) 

with A* — — £*, the energy of the excited KS orbital 
that the added electron joins. Eq{N) is the ground- 
state energy of the A'^-electron system, and E* {N + 
1) is the energy of the excited state of the {N + 
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FIG. 2: Exchange-correlation contribution to the excited 
electron affinity for a test system (See text). 

l)-electron system. The excited xc electron affinity 
accounts for relaxation when an electron is added to 
an A^-electron system, forming an {N -t- I)-electron 
excited state. (No correction is needed for ionization 
out of species I Q: the KS HOMO energy exactly 
equals minus the ionization potential.) 

Eqs. © and Q become exact in the weak inter- 
action limit where there is no coupling to other KS 
excitations. Dynamic correlation is completely miss- 
ing. Although derived for two one-electron atoms, 
the case for iV^ (2) electrons on fragment 1(2) (both 
neutral and open-shell) follows analogously. Eq. {T)) 
is an exchange approximation which may also be ob- 
tained from perturbation theory: One computes the 
Vac matrix element in an (A^2 + 1) electron Slater 
determinant, where the added electron is placed in 
an unrelaxed virtual orbital of the 7V2-electron KS 
potential. 

Before turning to TDDFT, we first give a sim- 
ple example to test Eq. 0. We consider one 
fermion in a one-dimensional parabolic well nT|, 

H = —^-^ + jx^- Excited states of two fermions in 
the well, interacting via a scaled delta-interaction, 
Vcc = \5(x — x'), may be numerically calculated 
at all interaction strengths A ^ . Figure HI com- 
pares the affinity calculated from energy differences 
(Eq. IHl, with Eq. 0: as expected, Eq. is exact 
for weak interaction. 

Returning to TDDFT, the diagonalization process 
above is effectively hidden in the structure of the xc 
kernel. Our approximation for /xc is motivated by 
the above analysis and also by the form of the inter- 
acting and KS density-density response functions. 

Because only single excitations appear in the KS 
response, in a SPA Xs has one pole at ojq = €2* — €2 
(<I>g of Eq. 01). The true density-density response 
however has two poles, one at the local excitation 
W2* , and the other at the CT excitation loct2' ■ Gen- 
erating an extra pole in Eq. and folding in the 
double excitation (0o,0o) ~^ (02*, ^o); requires a 
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dressed (i.e. frequency-dependent) SPA (see also 
Ref. 10), with the form [g|/Hxc(w)|g] = a+b/{uj—c). 
Note here that \q] = (j)o(r)(j)2* (r). This form is 
consistent with the earher diagonahzation analysis 
(yielding Eqs.lHl, which may be written 



{^g\H\%)-EHL- 



\{%\H\'^L 



-{{^d\H\<Pd) -Ehl) 



(9) 

A first approximation for the /xc matrix element 
then results directly from subtracting LUq = 62* — £2 
from the right-hand-side of Eq. ^ (c.f. Eq. |5J). Al- 
though this would give a huge improvement over any 
adiabatic approximation (see later), it lacks the cor- 
rect local Hartree-exchange-correlation effects. We 
now modify the approximation to incorporate these. 

As {<i>D\H\<pD) = ($,|i/|$9)-l-0(e-^) at large R, 
Eq. ||HJ) suggests a two-parameter model, 



2b|/Hxc('^)|g] ^a-ujq + 



(10) 



We set these parameters by requiring the solutions 
of Eq. with Eq. ^1 i.e. uj — a ± b, reproduce 
ATDDFT values for local xc effects and for lowest 
order polarization. That is, 



i {UJ2' + L0CT2' ) , 6 = i (W2- 



■ L0CT2' 



(11) 



where 



W2' 



LOg + A2. + D{1, 2*) - D{1, 2) , and 



WCT2- = ujg + AcT2^--+D{l+,{2-y)-D{l,2) 

K 

(12) 

(Note again ujq ~ £2* ~ £2 is the bare KS frequency). 
Here ATDDFT calculations on separated fragments 
determine all the quantities: (i) A2» is the AT- 
DDFT value for the excitation on fragment 2, mi- 
nus the KS frequency e2* — £2- For example, in 
a SPA, A2. - 2[q|/H + /4,ttI9]' ^^ere /^^^^ = 
'52£'xc/'5n|(r)(5n|(r'), with i?xc being a ground-state 
xc energy approximation. 

(I2 — /i) — ^ic*\ where we write 

- cj*(iV2 + 1) + e2-- Here A(N2) is 
the usual electron affinity of fragment 2, that may 
be obtained from ground-state energy differences be- 
tween the negative ion formed by adding an electron 
to fragment 2, and that of the neutral fragment 2. 
The frequency of the excited state of the (A^2 + 1)- 
electron ion 2, a;*(iV2 + l), is then given by ATDDFT 
performed on this negative ion. 
(iii) D(l,2) is the dipole-dipole energy between 
fragment 1 and fragment 2 in their ground states; 
Z?(l,2*) is that when fragment 2 is in the excited 



CT2* 



(ii) A, 

4c^ = A{N2) 



state; Z3(l+,(2^)*) is that between the positive 
donor 1 in its ground state and the excited acceptor 
state on fragment 2. The dipole moments can ei- 
ther be directly obtained from the ground-state DFT 
densities, or extracted from ATDDFT response on 
separated neutral or charged fragments. 

Thus we obtain the dressed kernel matrix element 
Eq. ijini) in terms of KS quantities and ATDDFT run 
on the separated neutral and charged fragments. 

As an example, consider two high lying excita- 
tions of the BeCl-|- molecule, that dissociates in its 
ground state to Be+ -I- CI. We consider the local 
excitation that dissociates to Be"*" (3s) -I- Cl(3p) and 
the CT excitation that is Be(3s) -I- Cl+(3p). We 
use the B3LYP functional for the ATDDFT pieces 
of the calculation, as programmed in the NWChem 
code j45j. Our approximate kernel yields the tail of 
the Be(3s) -I- Cl"'"(3p) CT dissociation curve to have 
the frequency (0.362 — l/i?)H, and the local Be"'"(3s) 
excitation to be 0.394H. These have the correct de- 
pendence up to OiXjB?) on the separation, and the 
asymptotes correspond to the generally reliable AT- 
DDFT values on the isolated Be-f and CI fragments. 
On the other hand, a purely adiabatic approxima- 
tion applied to the molecule yields only the local 
2s— » 3s excitation on Be+, with frequency 0.430H; 
giving only half the usual ATDDFT correction for 
a local excitation on an isolated fragment (see also 
Sec. lIV|) . and completely missing the CT one. Com- 
paring with the atomic spectroscopic data of NIST 
for the transition frequencies we find the local 
excitation to be 0.402H and the CT asymptote to be 
0.386H. 



IV. DISCUSSION AND OUTLOOK 

Eq. HlOfl. together with the parameters of Eq. H12|) . 
is our approximation for the kernel. The strong 
frequency-dependence is a consequence of static cor- 
relation in the KS system: our result does not go 
beyond the accuracy that ATDDFT has for usual 
local excitations, so it demonstrates the compli- 
cations static correlation creates for TDDFT. The 
static correlation is caused by the near-degeneracy 
of the HOMO and LUMO orbitals due to the step 
in the KS potential. Our result provides a practi- 
cal scheme to deal with this, constructed from AT- 
DDFT calculations run on the individual fragments, 
and builds in first-order molecular polarization ef- 
fects. We note that our result is only accurate for the 
tail end of the excited molecular dissocation curves 
(up to 0(l/i?^)). No xc effects across the long- 
range molecule are included. Molecular Feshbach 
resonances at frequencies higher than 1\ cannot be 
accurately described; these appear as KS shape res- 
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onanccs tunneling through the step. Limited in this 
way, and yet having comphcated form, our approx- 
imation highhghts challenges for a more complete 
description of dissociation. 

There is no exponential growth of the kernel ma- 
trix element with separation, that one would ex- 
pect from an adiabatic single pole analysis of long- 
range CT with localized (atomic)orbitals there 
the kernel would need to grow as in order to 
get a non- vanishing xc affinity. However, as shown 
here, for a long-range molecule consisting of open- 
shell fragments, this does not occur with the cor- 
rect delocalized molecnlsLT orbitals, provided our non- 
adiabatic approximation is used, as the CT state is 
"carried along" with the local excitation, through 
KS excitation out of the bonding orbital. 

We note that the lowest CT states in Ref. , and 
discussed in Sec.mihere. are a distinct case from the 
higher excitations considered here because there the 
underlying KS transition was nearly-degenerate with 
the KS ground-state itself: a single excitation to the 
antibonding transition (j}Q. By perturbative argu- 
ments, the KS response goes as the inverse frequency 
^ e^. The structure of the resulting kernel (Eq. Q 
in the present paper) is distinct from that for the 
higher local and CT excitations here (Eqs. ^| and 
I12|l . but also displays strong frequency-dependence. 

Any adiabatic approximation is stuck with un- 
physical half- electrons on each atom. The correc- 
tion to the bare KS energy would be half of what 
it should be for the local excitation (since the /xc 
matrix element involves the bonding rather than the 
localized orbital), and entirely misses the CT one. 

The step appears in the exact ground-state 



KS potential, and in orbital-dependent approxima- 
tions 13 . If instead LDA, or a GGA, was used for 
the ground-state potential, as in most calculations 
today, there is no step. The HOMO and LUMO 
are nevertheless also delocalized over the lo ng-r ange 
molecule. It is now well-established [s^, 113, 113 
that the dissociation limit of the molecule places 
(unphysical) fractional charges on each open-shell 
fragment in approximations such as LDA. The indi- 
vidual atomic KS potentials become distorted from 
their isolated form in such a way as to align the in- 
dividual homo's with each other, yielding molecu- 
lar orbitals that span the molecule, albeit unevenly. 
The HOMO and LUMO energies increasingly ap- 
proach each other as the molecule dissociates, yield- 
ing static correlation in the long-range molecule. 
Clearly, in light of the results here, an adiabatic 
kernel used in conjunction with such a ground-state 
potential, will fail for all excitations, and give un- 
physical fractional charges on the dissociating frag- 
ments in the excited states. Frequency-dependence 
is needed to fold in the double-excitation associated 
with excitation into the LUMO. Although the re- 
sult for the xc kernel of Ref. , and our results in 
Sec, mil are derived with the exact ground-state po- 
tential in mind, it may be interesting to explore the 
use of our kernel on top of an LDA ground-state po- 
tential, as much of the essential physics is captured. 
The results of this paper will be an important guide 
for likely modifications required in this case. 
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